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Abstract 

The anomalous dimension of the heavy-light quark current in HQET is calculated 
with three-loop accuracy, as well as the renormalized heavy-quark propagator. The 
NNL perturbative correction to /b//d is obtained. 
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1 Introduction 



Heavy quark physics progresses fast, both experimentally and theoretically. 
Extraction of fundamental Standard Model parameters from high-precision 
data from the 5-factories requires thorough theoretical understanding of 
strong-interaction effects. In particular, higher radiative and power corrections 
have to be calculated. 
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Understanding properties of hadrons with a heavy quark has been boosted by 
the Heavy Quark Effective Theory (HQET, see the textbook [f]). ft is an effec- 
tive fieid theory approximating QCD for probiems with a singie heavy quark 
when its 4- velocity v remains approximately constant, so that its momentum 
is p = mv + k, and the characteristic residual momentum of the heavy quark 
k is small compared to its mass m, and characteristic momenta of light quarks 
and gluons are also small. The HQET Lagrangian can be systematically con- 
structed as a series in 1/m. Operators of full QCD are also expanded in 1/m, 
coefficients are HQET operators with the appropriate quantum numbers. 

Here we are going to consider QCD heavy-light bilinear quark currents. They 
have numerous applications: vector and axial (with the anticommuting 75) 
currents describe heavy-to-light transitions with emission, tensor currents 
- the b ^ transition which appears at one loop in the Standard Model. 
The QCD heavy-hght current is jo — qoTQo — Zj{ai^^\fj,))j{fj), where F is 
a Dirac matrix, and Uf — ni + 1 is the full number of flavours including the 
heavy one. Its HQET expansion is 



where the HQET current is Jo = qo^hyo — Zj{as{ij))3{ij), — fh^, is the 
HQET heavy-quark field, and Oi are dimension-4 HQET operators having the 
appropriate quantum numbers. In order to have just one leading-order term 
in (1), we should treat f and J± — J^ — fv^ in T separately. 

The dependence of the matching coefficient C{fj,',ii) on the normalization 
scales /i' and /i is determined by the renormalization group: 



C{ii', 11) = C (m, m) exp 



7j(a^) das f Til^s) da^ 



2/3('^/)(aJ as J 2pM{as) a 



al (m) otl "(m) 



(2) 



where 



d log Z ^ as ^ f a 
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and a^J^^\fi) contains only rii light flavours. The HQET current anomalous 
dimension 7^ does not depend on the Dirac structure of the current F, and was 
known at one [2,3] and two [4,5] loops. The QCD current anomalous dimension 
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7j vanishes for the vector and axial currents (with the anticommuting 75) to aU 
orders. The QCD anomalous dimensions of the scalar and pseudoscalar (with 
the anticommuting 75) currents are known to be equal (modulo a sign) to the 
quark mass anomalous dimension, the latter is available at four loops [6,7]. 
The anomalous dimensions of the axial and pseudoscalar currents with the 
't Hooft-Veltman 75 are known at three loops [8]. Generic expressions for an 
arbitrary F are known at two [9] and three [10] loops. The matching coefficients 
C(m,m) = 1 + Cia^"'^(m)/(47r) + C2[a'^J^'\m) / (Att)]'^ + ■■■ contain no large 
logarithms; for all F, they are known at one [11] and two [9,12] loops. 

As one can see from (2), two- loop (next-to-next-to- leading order, NNL) correc- 
tions to C(m, m) should be used together with three- loop (NNL) anomalous 
dimensions and /5-function. Therefore, the results of [9,12] cannot be fully 
utilized without knowing 7^2- ^ The aim of this work is to calculate it. 

The 5-meson leptonic decay constant fs is one of the most important non- 
pcrturbativc parameters in B physics. It has not been measured experimen- 
tally; observation of the r^Uj. decay is a very difficult (though not 
hopeless) task for S-factories. Theoretical estimates vary widely. Improving 
their accuracy is extremely important. The HQET matrix clement <0|J(/i) |-B> 
at a low ^ mi, can be obtained by some non-perturbative technique, like 
lattice simulations (see recent review talks [13]) or HQET sum rules [14,15]. 
Direct lattice calculations with a heavy quark of mass ~ m;, will not be possi- 
ble in the foreseeble future, because they require the lattice spacing a <S i/rrih. 
A way to avoid the far extrapolation from m ~ rric (which can be simulated 
on present-day lattices) to m ~ is to use the lattice HQET, in which the 
heavy-quark mass m does not appear. Such simulations produce results nor- 
malized at a low scale /i ~ 1/a. They should be matched to the continuum 
HQET and run up to n — rrib using the three- loop anomalous dimension 7^, 
and then matched to QCD using the two-loop C{m,m) [9,12]. Therefore, cal- 
culation of 7j2 is one of the ingredients necessary to make lattice evaluations 
of fs more precise. Similarly, HQET sum rules produce results normalized at 
low 11 of the order of a typical Borel parameter, they also have to be run up to 
fj, — rrib. The knowledge of jj is also important for improving the accuracy of 
lattice and sum-rule calculations of the form factors of heavy-to-hght exclusive 
scmilcptonic decays B n, B ^ p (for soft tt, p) and exclusive rare decays 
caused by the b ^ s^y transition. 



^ 7j cancels in ratios of matrix elements, which were investigated in [9]. 
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2 Three-loop HQET propagator diagrams 



A systematic method of calculation of two-loop HQET propagator diagrams 
based on integration by parts [16] has been constructed in [5] . At three loops, 
this has been done in [17]. The algorithm has been implemented as a REDUCE 
package Grinder available at ht t p : / / www-ttp . physik . uni-karlsruhe .de/Progdata/ttpOO/ttpOO-01/ 

There are 10 generic topologies of three-loop propagator diagrams in HQET 
(Fig. 1). Grinder reduces all of them to linear combinations of 8 basis integrals 
shown in Fig. 2, with coefficients being rational functions of the space-time 
dimension d. 




f 9 




h i J 
Fig. 1. Topologies of three-loop propagator diagrams in HQET 

The first 5 basis integrals are trivial: 

j2 q2 

Bi = Ii , B2 = hh , -B3 = /a , -84 = I3— , i?5 = Is— i , (3) 

-'2 (-^2 

where 
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Fig. 2. Basis integrals 

r(i + 2n£)r"(i -£) 

(l-n(d-2))2n ' 
1 



r(l + n£)r"+i(l -e) 



" (n + 1 - nf )^ ((n + 1)| - 2n - l)^ r(l - (n + l)^) 
are the n-loop sunset HQET and massless integrals. 
The next two basis integrals are 

Be = 6-1/(1, 1, 1, 1, e), B, = hJ{l, 1, -1 + 26, 1, 1) . 



(4) 



We use recurrence relations for the two- loop integrals / and J (see [17]) to 
re-express them via /(1, 1, 1,1,1-1- e) and J(l, 1,2 + 2e, 1, 1). The integral 
7(1, 1, 1, 1, n) for any n was calculated in [18], and J(l, 1, n, rii, 712) for any n, 
rii, 712 - in [17] (Fig. 3). In particular. 



/(1,1,1,1,1 + e) 
X3F2I 



4r(l -e) 



r(l-h4e)r(l + 6£) 



9{d - 3){d - Ay [ r(i + 7£) 

3£, 3e, 6e 



\ 1 + 3£, 1 + 7£ 
J(l,l,2 + 2£,1,1) 



ij - r'(i + 3£)r(i - 3£) 
1 



r(i + 6£)r2(i-£) 



3{d - A){d - 5){d - 6){2d - 9) r(l + 2£) 
l,2-2£, l + 4£ ^\ 
V 3 + 2£,2 + 4£ / 

Such hypergeometric functions can be easily expanded in e in terms of multiple 



X 3^2 



(5) 
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( sums using the C++ library nestedsums [19]. Reducing them to the minimal 
basis [20], we obtain 



3^2 



, 1 + 3£, 1 + ?£ 



1 + 54C3£' - 513C4£^ + 54(25(5 + 28C2C3)£' 



+ 



3^2 



l,2-2£, l + 4£ 



1 ) = 2 + 6(-2C2 + 3)£ + 12(10(3 - IIC2 + 6)£' 

3 + 2£,2 + 4£ 

+ 24(-28C4 + 55C3 - 27C2 + 9)£3 
+ 48(94(5 - I6C2C3 - 154(4 + 135(3 - 45(2 + 12)£^ + 

Several more terms of these expansions can be found easily. 




(6) 



Fig. 3. Two-loop integrals / and J 

Finally, the most difficult basis integral can be found as following. We calcu- 
late the ladder diagram 1^(1, 1, 1, 1, 1, 1, 1, 1) (Fig. Ic, see [17]) using Grinder, 
and solve for B^. This diagram is convergent; at = 4 it is related [21] by 
inversion of Euclidean integration momenta to a known [22] on-shell diagram: 



1, 1, 1, 1, 1, 1, 1) = -5(5 + 12(2(3 + 0{e) . 



(7) 



A REDUCE program which expands in e results produced by Grinder can be 
obtained at [23]. 



3 Heavy-quark propagator 



Due to the non-abelian exponentiation theorem [24,25], the bare HQET heavy- 
quark propagator in the coordinate space can be written as 



S(t) = -ie(t) exp 



C 



9l 



it 



(47r)'^/2 V 2 



2e 



90 fit 



4£ 



{CaSfa + TpniSpi 



(47r)3<i/2 



{AttY V2, 

(it \ 
-j (clSpAA + CfTfUiSffi + CaTfUiSfm + (Tpn,)' Sfi 
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+ 



(8) 



Let's see how this happens exphcitely, up to three loops. 

Suppose we have calculated the one-loop correction to the HQET propaga- 
tor in the coordinate space (Fig. 4a). Let's multiply this correction by itself 
(Fig. 5). We get an integral in ti, ^2, t[, with < ti < t2 < t, < t[ < t'2 < t. 
Ordering of primed and non-primed integration times can be arbitrary. The 
integration area is subdivided into 6 regions, corresponding to 6 diagrams in 
Fig. 5. If the colour factors of the diagrams Fig. 4c, d were the same as that 
of the one-particle-reducible diagram (Fig. 46), i. e. equal to the square of the 
colour factor Cp of the one-loop diagram (Fig. 4a) (as in the abelian case), 
then the sum of the diagrams Fig. 46, c, d would be equal to | of the square of 
the one-loop correction Sp (Fig. 4a), as given by the square of the first term 
in the expansion of the exponent (8). In the non-abelian case, however, the 
colour factor of Fig. 4c differs from Cp by —CfCa/2, which is the colour factor 
of Fig. 4e. This is because when we reduce the colour factor of Fig. 4c to that 
of the reducible diagram Fig. 46, we get an extra term with the commutator 
[t", t''], which has the colour structure of Fig. 4e. Therefore, we should inchide 
the contribution of Fig. 4c with —CfCa/'2 instead of its full colour factor into 
the term SpA in the exponent (8). This part of the colour factor is called 
maximally non-abehan [24] or colour-connected [25]. Of course, the diagram 
Fig. 4e also contributes to Spa] the diagram Fig. 4/, with the one-loop gluon 
self-energy correction, contributes to Spi (quark loop) and SpA (gluon and 
ghost loops). 

^ ^ '^"^ 

a 6 c 

^ A\ 

d e f 

Fig. 4. One- and two-loop diagrams for the HQET propagator 

In the same way, we can consider the set of three-loop diagrams obtained by 
multiplying the corrections of Fig. 4a and /. We can imagine that this set 
is obtained from the one-particle-reducible diagram by allowing the gluon - 
heavy-quark vertices to slide along the heavy-quark line, crossing each other. 
These diagrams are said to contain two connected webs (or c-webs) [25]. Ev- 
erything is already accounted for by the product of the one-loop correction 
and the (Fig. 4/ part of) two-loop correction in the expansion of the expo- 
nent (8), except the contribution of Fig. 6a (and its mirror-symmetric), taken 
with the maximally non-abelian part of its colour factor. It contributes to 
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_ / \ _ t[ t'^ t 



Fig. 5. Exponentiation theorem 

the three-loop correction in the exponent. Similarly, out of all the diagrams 
with two connected webs Fig. 4a and e, only those of Fig. 66, c (plus their 
mirror-symmetric ones), and d contribute to Sfaa, with the maximally non- 
abclian part of their colour factors. This part appears, in the case of Fig. 6b, 
for example, when we commute matrices to obtain the colour structure of 
the reducible diagram; it is identical to the colour factor of Fig. 6h, equal to 






/ 
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Fig. 6. Three-loop diagrams contributing to the three-loop term in the exponent (8) 

Now we are going to consider all the diagrams with three gluons, both ends 
of which are attached to the heavy-quark line. They are said to contain three 
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c-webs, each of them is that of Fig. 4a. We decompose the colour factors 
of these 15 diagrams in the following way. We move the vertices along the 
heavy-quark line in such a way as to disentangle those c-webs. While doing 
so, we get extra terms from the commutators, having colour structures of 
the corresponding diagrams with the three-gluon vertex. These diagrams have 
fewer c-webs, which are more complicated. Finally, each colour factor can be 
expressed as a linear combination of three ones: Cp (3 c-webs of Fig. 4a); 
-CICa/2 (2 c-webs of Fig. 4a and e); CpCl/A (1 c-web of Fig. 6h). The first 
one occurs with the unit coefficient in all 15 colour factors. The sum of the 
corresponding contributions is just the term with the cube of the one-loop 
correction in the expansion of the exponent (8). The second colour structure 
occurs in the diagrams obtained by multiplying Fig. 4a and c; the sum of the 
corresponding contributions is contained in the product of the one-loop term 
and the two- loop one in the expansion of the exponent. We are left with the 
colour-connected parts of the colour factors (a single c-web contributions). 
They are present in the diagrams Fig. 6e, /, g (and its mirror-symmetric). 
They contribute to Sfaa in (8). 

We are left with the diagrams containing only a single connected web. Those 
of Fig. 6h and i have equal colour factors (this is evident if we close the 
quark line), they contribute to Sfaa- The colour factor of Fig. Qj is zero. This 
becomes clear if we close the quark line and write a thrcc-gluon vertex as the 
commutator (Fig. 7). The diagram with the four-gluon vertex (Fig. 6k) can be 
decomposed into three terms, with colour factors of Fig. 6h, i, j. The diagram 
Fig. 61 contains two-loop gluon self-energy corrections, including one-particle- 
reducible ones; it contributes to Sfaa, Sffi, Sfai, Sfu- The diagram Fig. 6m 
contains one-loop corrections to the three-gluon vertex, including one-particle- 
reduciblc ones (i. e., one-loop self-energy corrections to each gluon propagator 
in Fig. 4e); it contributes to Sfaa, Sfai- 



In order to obtain the heavy-quark field renormalization constant Zq, we 
should re-express this bare propagator in terms of the MS renormalized cou- 
pling as{ii) and gauge-fixing parameter a{ii) instead of Qq and qq:^ 



^ ao is defined so that the bare gluon propagator is [g/^u — (1 — cio)p^p,y/p'^] /p^. 




Fig. 7. Vanishing colour factor 



(47r)'^/2 




An 



ZaiasifJ,)) , ao = ZA{as{ii))a{ii) , 
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and require that the renormahzed propagator Sr{t; /i) — ZQ^S{t) is finite in the 
hmit e — > 0. After this re-expression, S{t) still has the exponential form with 
the same colour structures as in (8). Therefore, the renormalization constant 
can be written in the exponential form, too: 



Zq = exp Cf^Zf + Cf (CaZfa + TfUiZfi) (9) 

+ Cf {cIZfaa + CfTfUiZffi + CaTfUiZfai + {TFUif Zfu) + ■ ■ 

Hence the heavy-quark field anomalous dimension 7q = d\og Zg/dlog/j, has 
no Cp term at two loops and no Cp, C^Ca terms at three loops. 

We have calculated the HQET heavy-quark self-energy with the three- 
loop accTiracy, in an arbitrary covariant gauge. We have used GEFICOM [26] 
for automatic preparation of input data for Grinder. GEFICOM used QGRAF [27] 
for generation of diagrams. Then a specially written Mathematica program 
classified them into topologies, established the momentum routing required by 
Grinder for each topology. Dirac traces, index contractions and colour factors 
were performed with FORM [28]. Then Feynman integrals were calculated 
with the REDUCE package Grinder [17]. The whole process is completely 
automatic. We have checked that the bare propagator obtained from this self- 
energy, after Fourier transforming to the coordinate space, obeys exponentia- 
tion (8), exactly at any ao and d. The complete result for S{u;), as well as its 
expansion in e, can be obtained at [23]. 

Finally, we obtain the anomalous dimension 



^Q^2CF{a-3)-^ + CF 

4:71 



— + 4a 



179\ 32^ 
H Tpni 

6 y 3 

271 \ 



as 
47r. 



24a - f Cs 



23815\ 

216 ; 



+ 6CFTFnii-16C3 + 17) + CaTfUi (^-y a 



782 

96C3 + ^ 



as 



(10) 



It has the structure required by exponentiation. In the abelian theory without 
light fermions (Bloch-Nordsieck model), it equals the one-loop term exactly; 
it vanishes in the Yennie gauge a = 3 (which has been introduced for exactly 
this reason) . The two-loop term coincides with [5] . The three-loop term with 
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{TpUiY contained in the all-order large-n^ result [9]. The three-loop term 
has been obtained, by a completely different method, in [22]. Our calculation 
confirms this result. It is interesting to note that the difference of the HQET 
and QCD [29,30] quark field anomalous dimensions 7q — 7^ is gauge- invariant 
up to two loops [5] but not at three loops [22]. 

The renormalized heavy-quark propagator is 



1—1 Q. 



+ Cp (CpTpF + CaTfa + TfUiTFi) ( — 



+ Cf {CIT FFF + CfCa^ffA + CiT FAA 

V47r/ 



+ CfTfUi^ffi + CaTfUiEfai + {TlmfTFu) (|^)' + • • • (11) 



2 o , 

1537 „ 76 



Ej. = -4 , Tff = -2(20^ + 12(20 - 18C; 

(3\ 1537 76 

C2 + g) + (-C2 + 3) a - IOC2 - ^ , Tfi = 8C2 + y 

^FFF = -^Caa' + (24C3 + 8(2) + (-72(3 - 48(2) a + 72(3 + 72(2 - y , 

Effa = (4(3 - (2) + f-9(2 - l\ a" 



2 



/ 439 \ 

+ (^-124(3 + —(2 - 12j a + 264(3 - 315(2 + 



1537 



6 ' 

,4, , 29\ o / 1 , 19, 15, 197\ 2 
Efaa = -77(3 + C2 + — ]a^+[ -—U - -rCs + ^(2 + t^t « 



3^-^ 48/ V 16" 4' 4' 64 

•19^ 10 _ 175. 130. 187. 13921\ 



/19 10 175 130 187 13921\ 



a 



4603 2057 2167 4176353 
31(5 + 22(,(3 + ^(4 - ^(3 - -^(2 - 

128 . \ 179 



/ 128 \ 179 

Effi = (^32(3 - —(2] a - 24(4 - 200(3 + 120(2 + — , 

32, 2, 673\ 260, 2836, 



243 



_ 160 512 31232 
Efh - —(3 - -^(2 - 



For other normalization scales /i, it can be obtained by solving the renormal- 
ization group equation with (10). 

Analytically continuing the coordinate-space propagator from the half-axis 
t > to t — —IT, we have 
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+ CF{CAS^FA + TfniSU)[^y (12) 



+ Cf 

{^A^FAA + CpTpriiSppi + CaTfUiSpj^i + CFiTpniYSpn 

1537 __76 
24 ' ~ 3 ' 

29 



47r/ 



^F = 4, 5^^ = - a^ + (4C2-3)a-12C2 + ^, S^p, 



or 

'-^FAA — 



) 3 / 1 ^ 5^ 5^ 197\ 2 



48 

,19, 10,, 175, 304, 13921\ 



3 3 g „ '"9 576 / 

o-.. . 4603, 1089, 2533, 4176353 

+ 31C5 - 22C.C3 - + ^Ca - + , 

5];,, = 24C4 + 104C3 - ^ , ^^. = fC3 + ^, 

/^,.. 80. 673\ 556, 352, 184346 



4 Heavy-light quark current 



In order to calculate the anomalous dimension of the HQET heavy-light cur- 
rent, we consider the bare proper vertex function of this current. We put the 
light-quark momentum to zero. The heavy quark carries residual energy a; < 0; 
this is enough to ensure infrared convergence. The integrals are of the same 
propagator type as discussed in Sect. 2. We follow the same automatic proce- 
dure as in Sect. 3: the diagrams are generated by QGRAF (there are 237 of 
them), processed by the Mathematica program, traces and contractions are 
calculated with FORM, and integrals - with Grinder. The exact d-dimensional 
result and its ^-expansion are available at [23]. 

After re-expressing this bare vertex function via the renormalized quantities 
Ois{lj)-i it becomes times the renormahzed vertex, which is finite at 

£ — > 0. We have checked that the higher 1/e poles in Z-p are those dictated 
by the renormalization group. The renormalization constant Zj — Z^Zq'^ Z^^ 
(with Zq from Sect. 3 and Zq from [29]) is gauge- invariant, as it should be for 
a gauge-invariant current j. The resulting anomalous dimension reads 



I3 



-3^-4^ 



16C2 



+ Ca 4(2 



49 



10 



Tpui 



An) 



12 



C2 (^-80(4 - 36(3 + 64(2 - y ) 



176 



+ CfTfUi 



C3 + 



655\ 

9 

260 . 1451 \ 
448 470 \ 



^ _ ,152, 112, 512 
+ C,T,n,[—i,-—C,- — 



OisY 



(13) 



This is the main result of this paper. The two-loop term coincides with [4,5]. 
The three-loop term with [Tpniy is contained in the all-order large-n; re- 
sult [9]. All the rest is new. 

For the standard QCD {SUc{S) colour group) values Cp — 4/3, Ca — 3,Tf — 
1/2 eq. (13) becomes 



a. 



7j =■ 



-+ 



TT 



127 5 

1 ni 

72 36 



+ 



61 



+ 



43 



192 324 



ni + 



35 

1296' 



TT 

343 49 



89 83 475 
+ 72 ^=^+108''' ^^^"^^^ 

-— + (-3.04328 + 0.138889 ni ) (— 
(-12.941 + 1.55406 n; + 0.0270062 



TT 



TT 



(14) 



5 Application: Jb/Jd 



As a simple apphcation, let's consider the ratio fs/ fo- We have 

2 



1 



47r V 47r / 



F(=^)(m,) + of^ 
V rrir. 



where [9,12] 

ci = -2Cf , 



13 



CF(47r2log2-14C3-15C2 + 
+ Ca (-27r' log 2 + 5C3 + 5(2 - 



255 
~16 



871 
~A8 



+ Tpni (4C2 + ^)+Tf (-24C2 + ^) 



(15) 



and we have neglected corrections ~ {me/nib^ in c^'^ {mb^c are the on-shell 
masses of b, c). Similar formulas hold for Jb- Running of the HQET matrix 
element F^^\n) is governed by the anomalous dimension (13) with n/ = 4; 
F^^\mc) is related to F^^\mc) by the decoupling relation [12] 



89 



An 



+ 



(16) 



We obtain 




a>s [nib) ^/^o '^^ 

7,0 (iji Pi\ 56 4403 
""^ 2/5oUo /^oy' 75''' + 1875' 

8. 115 

r2o = Ac2 + ^2 = -C2 + ^ , 

8 21 o 152^ 254, 272392, 2071339 
= -7r log2 + -^C4-^C3 + ^^C2 + ^^^, 

/3-function and the anomalous dimension are for 4 flavours, Ac2 = — — 
CfTf (4C2 + fi), -22 = ICfTf (see (15), (16)). Numerically, 
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fn 




6/25 



l + 0.894(x-l)^£^ 



nib 

+ (3.788a;2 - 0.799a; - 2.448) 



TT 

.(4)r 



V ^ / V ^c,b , 



(1^ 



The coefficient of {as{mb)/TiY is about 5.5. Note tfiat tfie power corrections 
~ {rric/mbY ^-^e omitted from tlie NNL term, so tfiat tfiis formula does not 
apply at rUc ~ m^. Unfortunately, lattice and sum- rules estimates show that 
the Aqcb /iTic correction to fo is large. 



6 Conclusion 



We have calculated the HQET heavy-quark self-energy up to three loops, in an 
arbitrary covariant gauge. The exact rf-dimensional result and its e-expansion 
are available at [23]. Upon renormalizing the vertex function, we have obtained 
the heavy-quark field anomalous dimension (10), in agreement with [22], and 
the renormalized heavy-quark propagator (11), (12). We have calculated the 
vertex function of the HQET heavy-light current with zero momentum of the 
light quark, up to three loops, in an arbitrary covariant gauge. The results are 
available at the same site. From it, we have obtained the anomalous dimension 
of the heavy-hght current (13). This is the main result of this paper. We have 
derived the next-to- next-to-leading perturbative correction to fo/ fs (17). 
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